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Abstract 

We present the leading-logarithm (LL) final-state radiative effects for the exact O(a) 
YFS exponentiated (un) stable WW pair production at LEP2/NLC energies using Monte 
Carlo event generator methods. The corresponding event generator, version 1.12 of the 
program YFSWW3, wherein both Standard Model and anomalous triple gauge-boson 
couplings are allowed, generates 71(7) radiation both from the initial state and from the 
intermediate W + W~ state, and generates the LL final state W decay radiative effects. 
Sample Monte Carlo data are given for illustration. 
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The role of the final-state radiative (FSR) effects in the processes e + e~ — ► W + W~ + 
77,(7) ~~ > 4/ +71(7) at and beyond LEP2 energies is of considerable interest for the LEP2 and 
NLC physics programs 0, In this paper, we evaluate for the first time the possible 
interplay between the exact O(a) electroweak (EW) corrections and the leading- logarithm 
(LL) final-state radiative effects for these processes when the 72(7) radiation is realized 
according to the amplitude-based Monte Carlo event generator techniques described in 
Refs. [|], ||, wherein infrared singularities are cancelled to all orders in a by using the 
extension to spin 1 charged particles of the theory of Yennie, Frautschi and Suura for 
QED@. ' ' 

The final-state radiative effects are realized in the LL approximation using the calcu- 
lation of the program PHOTOS (Ref. 0) in which a non-radiative final-state process is 
used to generate up to two photons in the corresponding radiative process by iterating the 
structure function evolution equation^ for QED ||. The exact O(a) YFS exponentiated 
final-state W decay radiative effects will be published elsewhere IJ. In this connection, 
we note that we expect the non-leading 0(a) and higher order (0(a n ), n > 2) final-state 
radiative effects to be small, ~ 1% in the peak reduction effect for example, even for 
a "bare trigger" acceptance for the outgoing final charged particles. This has been found 
by the authors of Ref. 0, who analysed the effects of final-state radiation in Z decay 
in the naive exponentiated (exact and LL) 0(a) approximation and who estimated the 
corresponding size of the analogous effects in W decay, such as ~ 14% for the total peak 
reduction effect. Indeed, more recently, the authors of Ref. have made an independent 
cross check on their estimates of the FSR line-shape effects for e + e~ — > W + W~ —>■ 4f in 
Ref. 0. There they present an exact O(a) calculation of the process in the double-pole 
approximation (DPA), wherein they retain in the pole expansion [[K]] of the complete 
e + e~ — > 4/ amplitude, only the terms containing the double pole in the S-matrix at the 
complex mass squared, M 2 = Myy — iMw^w, where Mw,^w are the respective mass 
and width of the W boson, and where in the residues of the respective double poles they 
project the respective O(a) corrections to an appropriate on-shell point. Henceforth, we 
refer to the on-shell residue projected DPA as to the leading pole approximation (LPA), 
with more general applications in mind: for example, in a triply resonant process, the 
LPA would correspond to the triple pole terms in the respective S-matrix element with the 
residues projected to an appropriate on-shell point. In this gauge-invariant calculation, 
these authors find that the FSR peak reduction effect is ~ 14.4% for W + ^~' — > e + ^V e (^ e ), 
to be compared with their estimate of ~ 14% in Ref. . We will compare our results with 
those in Refs. |3| . We emphasise that our work differs from the work of Ref. in that we 
include the exact EW 0(a) corrections with YFS exponentiation in the production pro- 
cess and we actually calculate the effects of the FSR in the VF-pair production and decay 
process at LEP2/NLC energies; in Ref. 0, only the process v^v^ — ► ZZ — > e + e~ + v T v T 
is actually calculated, and a heuristic argument is used to estimate the corresponding 



1 To be precise an ansatz is provided, which reproduces the LL terms. It includes transverse degrees 
of freedom for the photon 4-momentum, assures coverage of the full phase space and rules of energy- 
momentum conservation. The photons' angular distribution is chosen to reproduce exactly the one of 
the soft photon limit. See Ref. [?| for more details. 
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results for final-state W-decay radiation. Thus, our calculations will also be a comment 
on the accuracy of these heuristic arguments in the presence of the YFS exponentiated 
exact 0{a) corrections to the H^-pair production process. Our work differs from that in 
Ref. H in that we include the YFS exponentiation of the exact O(a) production process 
in the Vy-pair intermediate state and the 0(a) 2 LL FSR whereas, in Ref. ||, the exact 
O(a) correction to the production and decay processes for the jy-pair in the leading 
pole approximation is calculated without exponentiation. The leading pole approxima- 
tion treatment of the attendant non-factorizable corrections in Refs. [11, 14 1 is also 



retained. The latter non-factorizable corrections have been shown pT| , [12| , |14| to be small 
and, as we explain in Ref. ||, when one works up to but not including C(f ^ 3L ) as we 
do, such effects may be dropped; which we do. Thus, although we start our calculation 
in Ref. || in the fermion-loop scheme | i~5| , when we focus on the 0(a) EW correction, 



we go to the leading pole part of the respective production amplitude. We also make the 
approximation of using on-shell residues for this double pole part, with which we then 
approximate the corresponding 0(a) EW correction. In our Monte Carlo event generator 
approach, we stress that the full off-shell phase space is always retained here. We improve 
our result by using the complete on-shell residues for EW corrections rather than their 
on-shell fermion-loop scheme representatives. Indeed, for the QED bremsstrahlung cor- 
rection we stress that, since the real photon has k 2 = 0, the corresponding running charge 
is the usual one. It can thus be shown that, in O(a), bremsstrahlung residues of the 
on-shell fermion-loop scheme are equivalent to those in the LPA; in both cases all infrared 
singularities are properly cancelled and not only the QED gauge invariance is preserved 
but also the full SU 2 l x U\ gauge invariance |4], ||. For this reason, in order O(a), in 
our final result, any reference to the fermion-loop scheme is purely pedagogical. What we 
arrive at is precisely the LPA, with full on-shell residues for the respective double pole 
approximation. Indeed, as the YFS expansion is not generally familiar, if one looks at 
Eqs. (1) and (2) in Ref. 0, which give the on-shell O(a) contributions to the YFS resid- 
uals Po and Pi, respectively, for the production process in the LPA, one may think that 
the lowest-order contribution to f3 , (3^ in the notation of Ref. 0, is not required either 
to be evaluated at the corresponding on-shell point as well. However, the right-hand side 
of Eq. (2) in Ref. 0], for example, involves the subtraction, from the corresponding on- 
shell O(a) bremsstrahlung cross section, of the product of the YFS real emission infrared 
function S || |J by the on-shell lowest order Born cross section; we need to stress that the 
YFS theory then forces the contribution to Po corresponding to this respective lowest-order 

Born cross section, /3^\ to be evaluated at the on-shell point as well. Thus, according to 
the YFS theory, Eqs. (1) and (2) in Ref. || are entirely equivalent to results in Refs. || 
for the production process, for the contributions up to and including terms 0(a). As can 
also be seen from the results in Refs. |§, these approximations are valid up to but not 
including jf-)- We then apply the YFS Monte Carlo methods of two of us (S.J. and 
B.F.L.W.) |l(|, as extended to spin 1 particles in Ref. J|, to arrive at the respective exact 
O (a) p ro d YFS exponentiated results realized in YFSWW3-l.il. Hence, we stress that, as 
far as the 0(a) correction to the production process under study is concerned, the results 
in Refs. P, O should be equivalent, in view of the many cross checks carried out by the 
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authors of Refs. [[[7], 18|, [19| on the two corresponding electroweak on-shell calculations 
used therein. 

More precisely, starting from the calculations in the program YFSWW3-l.il in Ref. |J, 
which feature the exact 0(a) pro d YFS exponentiated results for the process e + e~ — > 
W + W~ + n(j) — > 4/ + 72(7), we have interfaced the outgoing final state to the program 
PHOTOS j7|. The latter uses the structure function evolution equation for QED || to 
generate up to two final-state decay photons for each W according to the respective LL 
probabilities to radiate; here the corresponding angular distributions of the decay photons 
are all generated in accordance with this LL approximation as described in Ref. j7|. The 
net probability of the respective event is unchanged, i.e. the normalization of YFSWW3- 
1.11 is unaffected by this interface, which will be described in more detail elsewhere 0. 
We refer to the version of YFSWW3 that contains this final-state radiative interface to 
PHOTOS as YFSWW3-1.12 and it is available from the authors @. In what follows, 
we present some sample Monte Carlo data from YFSWW3-1.12 to look into the possible 
role of FSR in the presence of the 0(a) EW corrections. For definiteness, we focus here 
on the current LEP2 CMS energy of 190 GeV and on the SM couplings. The complete 
discussion of both LEP2 and NLC energies with the illustration of anomalous couplings 
will appear elsewhere ||. 

Specifically, in Figs. 1-8, we show the results obtained with YFSWW3-1.12 on the 
processes e + e~ — > W + W~ + 71(7) — > cs + iv^ t = e, /1, for the cosine of the W pro- 
duction angle distribution in the CM (LAB) system, for the W mass distribution, with 
both "bare" and "calorimetric" definitions of that mass, for the CMS lepton final energy 
distribution, for both calorimetric and bare definitions of that energy, and for the corre- 
sponding distributions of the cosine of the lepton decay angle in the W rest frame. We 
note the following properties of these results. First, concerning the W mass distributions 
in Figs. 1 and 5, we see that the respective average values of Mw are as given in Table 1. 
There, EW-ex denotes the exact 0(a) pro d calculation of EW corrections 0; EW-ap de- 
notes the approximate treatment of these EW corrections as given in Ref. |p0[ ; No EW 
denotes that the EW corrections other than the ones coming from LL (0(a 2 )) initial- 
state radiation are turned off. The calorimetric results are all closer to their respective 
NO FSR analogues than are the bare results, as expected. The effects of the FSR for the 
muon case are all respectively smaller than the corresponding results for the electron case, 
again as expected because of the smaller radiation probability for the muon. The size of 
the shift of (Mw) is generally consistent with the discussion in Ref. @, which deals with 
primarily the line shape (peak position and height); in detail we see that, in the presence 
of the FSR, at the level of our statistical errors, for an average quantity such as (Mw), all 
three calculations in the table are sufficient, as expected. With regard to the guesstimates 
made in Ref. concerning the peak reduction and the peak position shift, we see from 
the BARE curves in Fig. 1 that our result of 13.5% for the peak reduction in the e~ case 
(comparing the EW-ex curves with and without FSR) is in good agreement with the 14% 
guesstimate of Ref. f2j and with the 14.4% found in the recent 0(a) on-shell LPA results 
in Ref. ||. The ~ —57 MeV guesstimated in Ref. [[| for the corresponding peak position 
shift in the e~ case was recently updated to —77 MeV in Ref. M] for the [i case, the up- 
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79.725 ± 0.036 
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CALO 
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CALO 
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Table 1: The results of the 125 x 10 6 statistics samples (weighted events) (except for Born, where 
the sample is 540x 10 6 of such events) from YFSWW3-1.12 for the average value of Mw as computed 
with the levels of radiative corrections as indicated for both bare and calorimetric treatments of the 
final lepton. See the text for more details. 

dated expectation from Ref. || for the peak position shift is —39 MeV. For completeness, 
we note that the size of the peak reduction effect in the \i case has been found to be ~ 8% 
in Ref. whereas in Fig. 5 we find 7.6%, again showing good agreement between our 
results and those in Ref. 0. Indeed, to compare our results for the peak position shift 
with those just cited from Ref. ||, we have performed Breit-Wigner fits to our line shapes 
in Figs. 1 and 5 with the values, both fixed and floating, of the W width. The results of 
our fits are shown in Table 2. For comparison, the fits are done for two different mass 
intervals, from 78 GeV to 82 GeV, and from 76 GeV to 84 GeV, to illustrate the role of 
the wings of the resonance in the fits. From these results we find that the BARE peak 
position shifts are estimated using the narrow fit range as 80.168 — 80.240 = —72 MeV 
and 80.199 — 80.241 = —42 MeV for the e and \x cases, respectively. We also computed the 
shift in the average invariant mass (Mw) of the W in the narrow range from 78 GeV to 
82 GeV as another estimate of the peak position shift for the BARE trigger and we found 
— 81.5±1.4 MeV and — 43.9±0.9 MeV for the e and \i cases, respectively. Thus, both sets 
of estimators of the peak position shifts are in reasonable agreement with the results given 
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80.142 80.135/2.2547 


80.216 80.217/2.0831 
80.207 80.207/2.0892 


EW-ap 


78 - 82 
76 - 84 




80.166 80.168/2.2105 
80.142 80.134/2.2547 


80.216 80.217/2.0832 
80.207 80.207/2.0892 


EW-ex 


78-82 
76 - 84 


80.231 80.231/2.0443 
80.227 80.227/2.0372 


80.166 80.168/2.2105 
80.142 80.134/2.2547 


80.216 80.217/2.0832 
80.207 80.207/2.0892 


W ->■ /i 


Born 


78 - 82 
76-84 


80.241 80.241/2.0308 
80.250 80.250/2.0295 






No EW 


78 - 82 
76 - 84 


80.232 80.232/2.0342 
80.238 80.238/2.0307 


80.198 80.199/2.1196 
80.192 80.190/2.1481 


80.217 80.218/2.0731 
80.217 80.217/2.0845 


EW-ap 


78 - 82 
76 - 84 




80.198 80.199/2.1196 
80.192 80.190/2.1481 


80.217 80.218/2.0731 
80.217 80.217/2.0845 


EW-ex 


78-82 
76 - 84 


80.232 80.232/2.0343 
80.238 80.238/2.0307 


80.198 80.199/2.1196 
80.192 80.190/2.1481 


80.217 80.218/2.0731 
80.217 80.217/2.0845 



Table 2: The results of the Breit-Wigner line shape fit to the YFSWW3-1.12 MC sample for the 
W~ invariant mass distribution at Ecus = 190 GeV. The input values of the W mass and width 
were: M w = 80.23 GeV and T w = 2.03367033 GeV (this value was used in the IV-fix fit). The 
fits were performed for two W invariant mass M ranges - as indicated in the table. See the text for 
more details. 



in Ref. 00; in this connection, we recall the slight difference in beam energy between our 
studies (95 GeV) and those in Ref. |§ (92 GeV). Moreover, we see in Table 2 the same 
pattern of results as we see in Table 1: the FSR effects for the e case are more pronounced 
than those for the \i case; the calorimetric acceptance reduces the size of the FSR effects; 
the results are not very sensitive to the EW correction to the production process. If we 
compare the predictions with and without FSR for the EW-ex and no EW cases we get 
a measure of the modulation of the FSR on the EW correction. From the curves in our 
Figs. 1 and 5 and the respective plots of the Srad as defined in the figures we see that this 
modulation is as expected. Concerning the cosine of the W production angle distributions, 
we see the interplay of the exact EW corrections on the one hand and the FSR on the 
other. Further, we see that the approximate EW corrections of Ref. , while a definite 
improvement over the no EW corrections at all, are not sufficient to describe this interplay 
at the level of 0.5-1.0%. Similar remarks hold for the lepton energy distribution in the 
CM system, although the corresponding insufficiency is reduced to the level of ~ 0.3% 
for the BARE case, for example for electrons. Concerning the distributions of the cosine 



2 The fit mass shift and the peak position shift approach one another as the fit range approaches a 
zero size interval around the peak; a similar remark applies to the shift in the average mass relative to 
the range over which it is taken around the peak. 
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of the lepton decay angle in the W rest frame, we again see the importance of including 
both the EW corrections and the FSR in Figs. 4 and 8, for the electron and the muon, 
respectively. In all cases, the results for the muon, particularly the BARE results, are 
less affected by the FSR than are the corresponding results for the electron, as expected. 
We stress that our results in Figs. 1-8 are generally consistent with those in Ref. || keep- 
ing in mind that we treat the 0(a 2 ) LL FSR and the YFS exponentiated on-shell exact 
0(a) pro d production process, whereas Ref. treats only 0(a) corrections in our LPA, 
in which only on-shell residues are used. Indeed, in addition to the agreements already 
cited, we call the reader's attention to the normalization correction in Fig. 9 of Ref. : at 
the CMS energy of ^/s = 190 GeV, it is —11%, in very good agreement with our result in 
Ref. U, which is (1 + S proc i)(p w ) 2 — 1 = —11.1%, for the latter result, we have used Table 2 
in Ref. for the relative correction S pro d = —9.9% to the production process, and the 
result in Ref. |2TJ for the 0(a) correction to the leptonic partial width p w — 1 = —0.686%. 
In addition, we can note that, for the case of the rv T decay channel, our results are also 
consistent with those in Ref. || for the peak position shift and peak reduction effects. In 
view of our higher-order corrections, we find quite reasonable all the agreements noted 
here. A more detailed discussion of such comparisons will appear 0. We stress that we 
have arrived at our results through a MC event generator realization of our calculation, in 
which realistic, finite px, ri(j) radiation is incorporated in the production process on an 
event-by-event basis, whereas the results in Ref. are all semi-analytical. This enhances 
the significance of the general agreement of our results where they do overlap. 

The issue of whether the calorimetric results are more realistic than the bare ones 
appears to depend on whether one is talking about the muon or the electron^. For the 
electron, it is very difficult to separate the soft photons with energy < that are 
responsible for the FSR effects of the W line shape as discussed already in Refs. 0, ||; 
they are just a part of the electromagnetic calorimeter response in general, which is used 
to measure the electron energy. For the muon, the energy is usually measured by a muon 
chamber in which, in general, these soft photons are not present. Thus, for the electron, 
our calorimetric results are more realistic; for the muon, it is the other way around. In 
either case, we see that precision VT-pair production and decay studies need to take the 
interplay between the FSR and the EW corrections into account so as to obtain the most 
precise tests of the SM; our calculations in YFSWW3-1.12 offer an avenue to achieve that 
goal. 
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W — > e v e 




Figure 1: The invariant mass distributions of W~ reconstructed from its decay products, e~i7 e , 
four-momenta. In the left pictures the electron is treated exclusively ('bare' electron), while in the 
right pictures it is treated calorimetrically ('dressed' electron - its four-momentum is combined with 
four-momenta of all photons emitted within an angle of 5° around its direction). The input values 

are M w = 80.23 GeV, T w = 2.034 GeV. 
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Figure 2: The angular distributions of W~ reconstructed from its decay products, e~i7 e , four- 
momenta. In the left pictures the electron is treated exclusively ('bare' electron), while in the right 
pictures it is treated calorimetrically as defined in Fig. 1. 
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Figure 3: The distributions of the final state electron energy in the LAB frame. In the left pic- 
tures the electron is treated exclusively ('bare' electron), while in the right pictures it is treated 
calorimetrically as defined in Fig. 1. 
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Figure 4: The distributions of the electron decay angle's cosine in the W~ rest frame. In the left 
pictures the electron is treated exclusively ('bare' electron), while in the right pictures it is treated 
calorimetrically as defined in Fig. 1. 
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Figure 5: The invariant mass distributions of W~ reconstructed from its decay products, n^v^, four- 
momenta. In the left pictures the muon is treated exclusively ('bare' muon), while in the right pictures 
it is treated calorimetrically as defined in Fig. 1. The input values are Mw = 80.23 GeV, F\y = 
2.034 GeV. 
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Figure 6: The angular distributions of W~ reconstructed from its decay products, /U - !^, four- 
momenta. In the left pictures the muon is treated exclusively ('bare' muon), while in the right 
pictures it is treated calorimetrically as defined in Fig. 1. 
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Figure 7: The distributions of the final state muon energy in the LAB frame. In the left pictures 
the muon is treated exclusively ('bare' muon), while in the right pictures it is treated calorimetrically 
as defined in Fig. 1. 



15 



W — > fi 

"BARE" "CALO 




Figure 8: The distributions of the muon decay angle's cosine in the W~ rest frame. In the 
left pictures the muon is treated exclusively ('bare' muon), while in the right pictures it is treated 
calorimetrically as defined in Fig. 1. 
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